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Abstract 

This report describes a cantilever controller for magnetic resonance force microscopy (MRFM) based on a 
field programmable gate array (FPGA), along with the hardware and software used to integrate the controller 
into an experiment. The controller is assembled from a low-cost commercially available software defined radio 
(SDR) device and libraries of open-source software. The controller includes a digital filter comprising two 
cascaded second-order sections ("biquads"), which together can implement transfer functions for optimal 
cantilever controllers. An appendix in this report shows how to calculate filter coefficients for an optimal 
controller from measured cantilever characteristics. The controller also includes an input multiplexer and 
adder used in calibration protocols. Filter coeSicients and multiplexer settings can be set and adjusted by 
control software while an experiment is running. The input is sampled at 64 MHz; the sampling frequency 
in the filters can be divided down under software control to achieve a good match with filter characterisics. 
Data reported here were sampled at 500 kHz, chosen for acoustic cantilevers with resonant frequencies near 
8 kHz. Inputs are digitized with 12 bits resolution, outputs with 14 bits. The experiment software is 
organized as a client and server to make it easy to adapt the controller to different experiments. The server 
encapusulates the details of controller hardware organization, connection technology, filter architecture, and 
number representation. The same server could be used in any experiment, while a different client encodes 
the particulars of each experiment. 
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FIG. 1: MRFM experiment 



I. INTRODUCTION 

This report describes a cantilever controller for magnetic resonance force microscopy (MRFM) 
based on a field programmable gate array (FPGA). The controller is assembled from a low-cost 
commercially available software defined radio (SDR) device, libraries of open-source software, and 
some additional software that we have written. In addition to the controller itself, we also describe 
software that integrates the controller into an experiment. We show how to calculate filter coeffi- 
cients for an optimal controller from measured cantilever characteristics. Finally, we present data 
comparing calculated and measured controller performance. 

Force microscope cantilevers are mechanical oscillators that require feedback control. In MRFM, 
cantilevers sense the magnetic force exerted by spins in the sample [H [2] . Cantilevers for detecting 
such small forces must have low spring constants and high quality factors, but these result in large 
deflections, long ring-down times, and narrow bandwith, which hinder image acquisition. These 
effects can be mitigated by feedback control without affecting the signal-to-noise ratio for force 
measurement [3] . In the apparatus shown in Fig. [T| the cantilever controller would be connected 
between RX and TX in the lower right corner (RX is the controller input). 

Cantilever controller design is discussed in [3] IH El [6]; details of our design appear in the 
appendix. Early controllers used analog circuits including operational amplifiers jH [5]; recent 
controllers use specialized microprocessors called Digital Signal Processors (DSPs) O [7] . In this 
report we describe a new digital cantilever controller based on an FPGA. 

An FPGA contains an array of identical digital logic and memory elements which can be con- 
nected together under software control to form, in effect, a custom integrated circuit. FPGA-based 
signal processors can have several advantages over microprocessors (including DSPs). Compu- 
tations can be highly parallel (instead of executing a sequential instruction stream); functions 
executing in parallel do not affect each other's speed, so there can be greater throughput even at 
lower clock speeds. Latencies can be lower because each input and output can be connected directly 
to the FPGA (instead of connecting them all through a shared bus). Precision (word length) can 
be tailored to computational requirements (instead of relying on a few fixed word sizes). 
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FIG. 2: System hardware and software 



This controller accommodates input and output signals up to 2 volts peak-to-peak (-1V to -|-1V), 
digitized with a resolution of 12 bits (input) and 14 bits (output). The sampling frequency can be 
64 MHz, or can be decimated (divided down) under software control. The controller demonstration 
described here samples at 500 kHz, chosen for acoustic cantilevers with resonant frequencies near 
8 kHz. The controller includes an infinite impulse response (HR) digital filter comprising two 
cascaded second-order sections (biquadratic sections, commonly called "biquads"). The signal and 
the filter coefficients are represented by 24-bit integers. The experiment software translates floating 
point to appropriately scaled 24-bit integers. 

In addition to the filter, the controller also includes an input multiplexer and adder that selects 
and optionally adds two input signals in all combinations. This is helpful for some calibration 
protocols. 

The controller characteristics (filter coefficients, multiplexer settings etc.) can be set and ad- 
justed by control software while an experiment is running. 

We expect that this controller will acccommodate future microscopes in our laboratory with 
modest effort, and we believe it could be used elsewhere as well. 



II. HARDWARE 

Figure [2] shows the cantilever controller system, including the two main hardware components 
(solid boxes): the host computer and the USRP (Universal Software Radio Peripheral), which acts 
as the controller. The USRP and host communicate over a USB 2.0 (Universal Serial Bus). The 
USRP is in its own small case (21 cm wide, 5 cm high, 14 cm deep) and the bus runs through a 
cable (up to 5 m, longer with extensions), so it is easy to connect the USRP to a different host. 



A. USRP (Controller) 

All of the controller hardware including the FPGA, data converters, timing, signal conditioning, 
and connectors are included in the USRP [H |9]. The USRP is designed to support software- 
defined radio, in particular the open-source GNU Radio project [10 . Its design is open source; 
its schematics are included in the GNU Radio software distribution. The USRP is commercially 
available. 

Figure [3] shows the signal path through the USRP: the signal travels from left to right (horizontal 
solid lines), digitized control parameters from the host computer enter on the USB bus (vertical 
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dotted lines) and the digitized acquired signals for the host computer exit on the USB (vertical 
solid lines). The multiplexer, decimator, and filters execute on the FPGA, an Altera Cyclone 
EP1C12, which provides 12,060 logic elements (each comprising several logic gates) and 239,616 
bits of storage [TT] . The data converters on the USRP are included in two Analog Devices AD9862 
codecs, each with two analog to digital converters (ADC) and two digital to analog converters 
(DAC) [12]. The ADCs have 12 bit resolution; the DACs have 14 bit resolution. The AD9862 also 
includes programmable gain amplifers (PGA) at the input of each ADC and the output of each 
DAC. The two blocks labelled "Codec in" and "Codec out" in Fig. [3] represent different parts of a 
single AD9862 device (there is one additional DAC with PGA which we do not use). 

The ADC output pins are connected directly to FPGA input pins, so the digitized signal reaches 
the FPGA with less latency than would be possible if the ADCs were connected through a bus 
shared by several devices (as is common in microprocessors, including DSPs). 

The ADCs sample at 64 MHz and the FPGA is clocked at 64 MHz. The DAC samples at 128 
MHz. The FPGA program includes a decimator {^N in Fig.js]) which can divide down the sampling 
frequency used by the filters. The precision required for the filter coefficients increases as the 
sampling frequency grows relative to the resonant frequency of the filter. A programmable sampling 
frequency makes it possible implement filters that work well over a wide range of signal frequencies, 
even when the filter coefficients are represented with limited precision. A good compromise for a 
filter with resonance near 8 kHz divides the 64 MHz clock by 128 to sample at 500 kHz. 

The FPGA, codecs, and USB hardware are on the USRP motherboard. The USRP also in- 
cludes two daughterboards, LFRX and LFTX (Low- Frequency Receiver and Transmitter, where 
low frequency means below 30 MHz). These provide connectors and AD813x differential amplifiers 
at the inputs and outputs for buffering and signal conditioning [13_. The amplifier circuits are 
fixed at unity gain and 50 ohm impedance. The daughterboards have no programmable elements 
(daughterboard components do not appear in Fig. |3]). 

The signal path shown in Fig. |3]only uses half the USRP: one of two AD9862 codecs, and two 
of four daughterboards. Another independent signal path could be accommodated in the same 
USRP. Both signal paths would execute in parallel so both would run at full speed. 



B. Host computer 

The host computer stores the controller program that runs on the FPGA. Before an experiment, 
a program on the host computer loads the controller program into the FPGA. During an experi- 
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ment, the host computer executes experiment control software that can set and adjust controller 
parameters and acquire the measured signals. 



III. SOFTWARE 

A. Controller software 

The controller program comprising the multiplexer and digital filter executes on the FPGA in 
the USRP. 

The multiplexer selects among input signals to connect to the digital filter. It can connect the 
cantilever position signal or a reference signal (RX or Ref in Fig. [3|, or neither, or their sum. All 
four combinations are used in our calibration protocols. 

The filter consists of two second-order stages, also called biquadratic stages, commonly called 
"biquads". Each stage implements a transfer function represented by a ratio of second-order 
polynomials, with three coefficients in both the numerator and denominator. The signals and 
coefficients are 24-bit integers. Using biquad stages requires less precision than a single higher- 
order filter. 

Two biquad filter stages are sufficient to implement transfer functions for optimal cantilever 
controllers derived by the method of the appendix. The transfer functions are computed and then 



translated to biquad stages by software on the host computer (section IIIB). 

In addition to the controller itself, the FPGA also executes code that communicates with the 
host computer over the USB. The bus control code, and several modules used by the controller, 
are from the open-source GNU Radio project |10j . 

The controller program is written in the Verilog hardware description language (HDL) [14j. The 
several Verilog source files are compiled or synthesized into a single binary file called a hitstream 
which is stored on the host computer and loaded into the FPGA on command. An FPGA can be 
reprogrammed by loading another bitstream. The host computer can store several bitstreams that 
provide different controller capabilities and load them into the FPGA under program control. 

The word length is the number of bits used to represent numbers in a program. An FPGA 
program can specify any word length at each stage in each computation, unlike a conventional pro- 
cessor or DSP where only a small fixed set of word lengths is supported by the machine arithmetic. 
FPGA capacity is limited, so an important design task is to choose the minimum word lengths 
that provide sufficient accuracy. Different word lengths can be chosen for the signal path (solid 
lines in Fig. |3]) and for the filter coefficients (dotted lines near the center of Fig. |3]). 

If the word length of the signal path is too small, the filter exhibits a nonlinear input /output 
relation with a threshold at small amplitudes and saturation at large amplitudes. With a 16-bit 
signal path, the linear range is only from about 0.1 V to 0.5 V input amplitude, which is too little 
dynamic range for our experiments. The word length required for the filter coefficients increases 
as the sampling frequency grows relative to the resonant frequency of the filter. If the word length 
of the coefficients is too small, the actual resonant frequency of the filter differs from the intended 
(computed) value. A good compromise for a filter with resonance near 8 kHz divides the 64 MHz 
clock by 128 to sample at 500 kHz. Under these conditions, 16 bit coefficients can result in an 
actual resonant frequency that differs about 16 Hz from the intended value, which is too much for 
our experiments. We found that 24 bits were adequate for both signal path and coefficients; some 
intermediate results use word lengths up to 50 bits. We found the acceptable word lengths by trial 
and error, although did do some simulations in MATLAB [TB] to confirm that the symptoms we 
observed were explained by the limited word length, rather than programming errors or some other 
cause. 
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We have only very limited ability to observe and debug synthesized Verilog code executing on 
the FPGA. Therefore we do most debugging in simulation. First we check the design in MATLAB, 
as described above. Then on the host computer we write the Verilog code and use the Icarus 
Verilog compiler to build an executable simulation (TB]. We execute the simulation on the host and 
use the GTKWave viewer on the host to show its results in a display that resembles an oscilloscope 
or logic analyzer [T^ . Finally we use a synthesizer to compile the debugged Verilog to a bitstream 
that can be executed on the FPGA |18j . 

The bitstream for this controller program uses 2,765 logic elements and 66,304 memory bits in 
the FPGA, about one fourth of the available resources. Another independent signal path could 
be accommodated in the same FPGA, or higher-order transfer functions could be implemented by 
adding more biquad filter stages. 

The controller software is open source and is available from the project web site |19j . 



The measured performance of the controller is discussed in section IV 



B. Host computer software 

The experiment software including calibration protocols, data acquisition and storage, calcula- 
tion of the controller transfer function, etc., executes on the host computer. 

The host software is organized as a client and server. The client comprises most of the exper- 
iment control, while the server is a small component that copes with the details of the controller 
and its interface. For example, the client calculates the transfer function for the optimal can- 
tilever controller as described in the appendix. The transfer function is represented as a ratio of 
polynomials with four numerator and four denominator coefficients. The server transforms this 
transfer function for two cascaded biquad sections by factoring the polynomials, sorting and col- 
lecting factors, then multiplying them in two groups. It also replaces floating point coefficients 
with scaled 24-bit integers. It loads the integer coefficients into particular registers in the FPGA 
using the USB. The client software does not depend on controller details such as filter architecture, 
hardware organization, number representation, and connection technology; these only affect the 
server. Therefore the same server could be used with different experiments; the particulars of each 
experiment are encoded in the client. Thanks to this organization, our client software is largely 
unchanged from earlier versions that used completely different controller technology based on a 
DSP. 

Client and server communicate by sending text messages over a TCP/IP socket connection. In 
Fig. [sj both execute on the same computer. They can also execute on different computers (we have 
used both configurations). 

The server also performs data acquisition, acquiring the stream of digitized samples from the 
FPGA at two points in the signal path: the outputs of the multiplexer and the digital filter (labeled 
"sum" and "filter" in Fig. [s]). The server can optionally display these data as a spectrum or time 
series on an oscilloscope-like display, as shown in Fig. |4| 

The server is written in Python and uses libraries from the open source GNU Radio and SciPy 
projects [To] [201121]. 

Our usual client is a large experiment control program written in Lab View [22^ . We also wrote 
a small client in Python for diagnostics and testing. Both clients can run at the same time, because 
both open and close the socket connection each time they send a message. 

The host software is open source and is available from the project web site |19] . 
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FIG. 4: Server program oscilloscope-like display 



IV. PERFORMANCE 

We compared a computed transfer function to the transfer function infered from measurements 
of the analog signals at the input and output of the controller, and found close agreement (Figs. [S] 
and [6]) . For this test we used a controller transfer function that we had computed for a typical 
MRFM experiment, by the method of the appendix, where the cantilever had resonant frequency 
near 8000 Hz and an uncontrolled resonant quality Q near 10,000. The transfer function is described 
by the four b (numerator) and four a (denominator) coefficients in table ^ 





7.026189 • 10-5 


1.027999 • 10-4 


-5.927540 • 10-5 


-9.181339-10-5 




1.000000 • 10° 


-2.848528 • 10° 


2.708790 • 10° 


-8.588522 • 10-^ 



TABLE I: Controller transfer function (index i ranges from to 3, left to right) 

Plots of the computed transfer function appear as solid lines in Figs. [5] (amplitude) and [6] 
(phase). At the resonant frequency, the ratio of the output amplitude to the input amplitude is 
about 1.4, and the output phase lags the input phase by about 80 degrees. The server software 
factors and scales this transfer function into two cascaded biquad sections with 24-bit integer 
coefficients (table 
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Section 1 
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2293 


-32865 


35158 


-49146 







-4194304 


8339278 


-4187298 


-4194304 


3608314 






TABLE II: Controller filter coefiicients (in each section index i ranges from to 2, left to right) 
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We measured the controller transfer function with a lock-in amp lifer [23]. The SINE OUT signal 
from the lock-in was applied to the input to the controller (RX in figs. |2]and[3]). The experiment 
software commanded the lock-in amplifier to generate sine waves from 7700 to 8300 Hz at 100 Hz 
intervals. The sine wave amplitude was adjusted to 0.100 V. The output from the controller (TX 
in the figures) was applied to the SIGNAL IN A-I input of the lock-in. At each frequency, the 
in-phase and quadrature components X and Y of the controller output were calculated by the lock- 
in amplifier, with a time constant of 3 sec, and the experiment software calculated the amplitude 
and phase of the transfer function from X and Y samples. The measured transfer function points 
appear as crosses in Figs. |5] (amplitude) and [6] (phase) . The errors (differences between calculated 
and measured values) appear in the plots below the transfer functions. The amplitude ratio error 
is less than 0.002, and the phase error is less than one degree, over the entire measured frequency 
range. These errors are acceptable in our experiment. 

The computed and measured transfer functions have a phase lag of about 80 degrees at the 
resonant frequency. The computed transfer function includes phase compensation to account for 
the phase lag contributed by various components in the entire signal path (appendix. Fig. [9] and 
section 3e|. The lag contributed by the fixed latency of the digital controller was estimated by 
applying a 50 Hz square wave to the input and observing (on an oscilloscope) the delay until the 
output starts ringing. This showed a latency of about 2 usee, similar to the sampling interval (at 
the 500 kHz sampling frequency), so the digital controller contributes about 6 degrees of lag at 
8000 Hz. 

The gains of the PGAs (Fig. [s]) were set to to 1, nominally. This results in (nearly) equal 
amplitudes for RX and TX, when a controller with unity gain is loaded. The PGAs are not 
precisely calibrated so the overall gain is not exactly 1.0. It is possible to trim the output PGA (via 
software) to obtain an overall gain very near 1.0, but we choose instead to simply keep track of the 
observed overall gain in our experiment software. The observed overall gain in these measurements 
was 1.038. 



V. APPLICATION 

We used the controller in an experiment with a cantilever with a resonant frquency of 8347 Hz 
and an uncontrolled Q of about 10,000. Figure [7] shows the one-sided power spectrum of the noise 
measured (by a spectrum analyzer [2^) from the controlled cantilever with several different con- 
troller transfer functions, with controlled Q values 75, 100, 200, 300 and 400. These measurements 
show the expected results. 



VI. DISCUSSION 

We anticipate adapting the platform described here to different control schemes. The following 
paragraphs describe one that we are considering. 

A goal of MRFM is to achieve atomic-scale resolution by resolving individual nuclear spins. A 
promising approach to this goal would use cantilevers that oscillate at radio frequencies (1-10 
MHz) p5] . Such high cantilever frequencies press the performance limits of available digital control 
technologies. An alternative uses heterodyne control, where the central controller elements process 
signals at much lower frequencies than the cantilever oscillation |26]. It confers another potential 
advantage because a heterodyne controller might also act as the lock-in detector for the MRFM 
signal, eliminating the separate lock-in amplifier we now use in our experiments [7]. 

In earlier work, we demonstrated a digital heterodyne controller with an emulated radio fre- 
quency cantilever |26j (our emulation was a real-time digital simulation with analog inputs and 
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outputs that resembled the signals expected from a radio frequency cantilever). We have begun 
experimenting with an FPGA implementation of a heterodyne controller [27. . The hardware and 
host computer software is the same as described here. The only difference is the FPGA code for 
the controller, which implements a different block diagram than the one shown here in Fig. [3j 
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APPENDIX: CONTROLLER TRANSFER FUNCTION 

Here we show how to calculate filter coefficients for an optimal controller from measured can- 
tilever characteristics. We derive the controller transfer function, following the theory of optimal 
feedback control of force cantilevers in |^ . After reviewing briefly the continuous-time specification 
of that controller, we derive a discrete-time equivalent. Practical hardware and software issues, 
along with compensation of parasitic effects, are treated. 



1. Continuous-time Control 

In Fig. [8j a linear dynamic feedback control element is added to the open loop system. The 
dynamic relationship between the forces u applied to the cantilever tip and the motion x of the tip 
is represented by the transfer function G{s). In addition, / is the force signal to be detected, w is 
additive process noise force, v is measurement noise added to the actual position of the cantilever, 
and H(s) is the transfer function of the feedback controller. 
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Many techniques are available for the design of H{s). All require knowledge of the physical 
characteristics of the cantilever, the noise processes w and v, and a design goal for controller 
performance. Below, the techniques of optimal control and estimation are combined to determine 
the continuous-time transfer function of the controller Hoc{s) which relates measured motion of 
the cantilever y to the control component of the input force u (Eq. (31) of [3]). 



where 



= -kaf3{a + (3) + —— , 



a/3 _ 9 I I _2_ 



a + /3 + ^ 



2 Q 



oc 
,2 



Qo 



Q 



where k, uJn, and Q are the stiffness, natural frequency, and resonant quality of a single-mode 
cantilever model. The parameters a and /? are determined from the control and estimation re- 
quirements. 

The value of a is defined by the need to keep the system stationary within acceptable deviations 
from X = 0, while using acceptable amounts of control u. 



a is a function of the maximum allowed variances of the displacement X = E[xf^]max and the 
control force U = E[v?\rnax- The definitions of ^ and a indicate that the relevant physical parameter 
is the ratio of the maximum control force to the force that would be caused in a spring subjected 
to the maximum displacement. 

The estimator parameter (5 represents the quality of the data used to estimate the system state, 
and is determined by the relative amounts of process and measurement noise. 



where constants W and V are the power spectral densities of stationary, zero-mean, Gauss-Markov 
random processes. Large measurement noise means that the observations are weighted less in the 
estimate of the state, while small measurement noise allows large gains in the estimator. 
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The controller has two complex poles and a real zero. As the allowable control effort is increased 
(greater a), or the measurement noise is reduced (greater (5), the natural frequency of the controller 
rises and its Qoc is diminished. Again, in the limit as (a + /3) approaches infinity, Qoc approaches 

y/2/2. 

Ultimately, the control accuracy can be no better than the inferred state. The control value, 
(u^), approaches infinity as a — > oo. However, even when large control effort is available, little is 
gained by allowing a to be arbitrarily large. The poles of the closed loop system are the union 
of the poles of the optimal controller and estimator. Since poles of smaller magnitude will tend 
to dominate the dynamic response, a practical guide is to make a smaller than /?, but of similar 
magnitude. 

While the closed-loop system is of fourth order, it is often useful to describe the controller design 
goal in terms of the effective resonant quality Qd of the dominant closed-loop poles. Assuming 
that a is smaller than /3, the value of a corresponding to the closed-loop quality Qd is given by 
(Eq. (16) of a) 



a 



1 

Q 



/2Q2 
2Q2 



d 



1 



1 



1 

Q' 



(A.4) 



Alternately the design goal can be stated in terms of the closely related quantities: the bandwidth 
Aujd = uJn/Qd, or the damping time Td = 2Qd/oJn- 

Finally, the performance of the optimally controlled closed loop system in terms of the mean- 
square cantilever displacement and control force can be predicted from (Eqs. 38-39 of [3]) 



Vf3[ 



+ ^-^ + aP + 2]['-^ + M 



(A.5) 



and (Eq. (39) of [3]) 



(n2) _ V{kaPf 

^ ~ 8(a + ^)(a2 + ?f +2) 



+ 2a^(5 + 2a^ + 20^ + -a^0^ + 8-h 



Q 

8a + 16P 8 
Q3 +Q4 



(A.6) 



Since /? is a function of V, it is useful to express V in Eqs. (|A.5|) and (A.6) in terms of the 

1 



(experimentally measured) uncontrolled (open loop) rms amplitude {y^ 



oil 



V 



(/34 + 4/32)Q2 + 4(4^ +2p + /33)Q + 4/32 



(A.7) 



where AuJv is the effective measurement noise bandwidth. 
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FIG. 9: Implementation block diagram 



Discrete-time Control 



The digital controller shown in Fig. |3] implements the feedback element H{s). The digital 
controller periodically samples the voltage signal corresponding to the measured cantilever motion, 
computes the appropriate control law, and imposes the result through a zero-order hold at its 
output. In our case, the control law is a discrete-time equivalent of the continuous-time transfer 



function in Eq. (A.l). 

In general, the control law implemented by the digital controller can be any linear causal time- 
invariant difference equation of the form 

N M , 

where y{n) and x{n) are samples of the output and input respectively, and at and bk are constants. 

Before we derive the difference equation several hardware and software issues must be consid- 
ered. 



3. Implementation Considerations 

Figure |9] shows the complete cantilever control system, including the elements of cantilever 
motion sensing, force actuation, conversion and amplification. Below, we consider the specifications 
and characteristics of each block. Throughout, we assume that the A/D and D/A converters are 
linear with inversely related gains: Kad = ^/K^a- 



a. Sampling Rate 

Faithful discrete-time reproduction of the continuous-time controller requires that the sampling 
rate be large in comparison with the natural frequency of the system. A good rule-of-thumb is to 
select a sampling frequency at least 30 times larger than the cantilever natural frequency. Typically, 
our sampling rate is 0.5 MHz for a natural frequency of 7500 Hz. 



h. The interferometer 



The fiber optic interferometer converts the displacement of the cantilever into a proportional 
voltage signal. A typical sensitivity for our interferometer is Ki = IQ x 10^ V/m. The interferom- 
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eter bandwidth is limited to approximately 100 kHz. 



c. Force Actuator 

Magnetic actuation provides the restoring feedback force to the cantilever. The force is propor- 
tional to the current passing through the control coil. The actuator can be calibrated by driving the 
coil (open loop) with a sinusoidal current of known amplitude at the cantilever natural frequency, 
and measuring the resulting cantilever motion. From that measurement and the known stiffness, 
quality, and natural frequency of the cantilever, the actuator sensitivity Kc is determined. 

A typical actuator gain is i^^c = 35 x 10^^^ N/A. 



d. Amplification 

Two stages of amplification are used to match the dynamic ranges of the interferometer and 
coil actuator to the input/output ranges of the digitial controller. 

Suppose that the full-scale ranges of both the A/D and D/A converters are Vfs- A simple 
strategy is to select the interferometer pre-amplifier gain Kpa such that the predicted rms value of 
the controller input is one third of the A/D converter range: 



Kpa = , ^ ,^ (A.9) 

Then select the current amplifier gain Ka such that the predicted rms value of the controller 
output is one third of the D/A converter range: 

(A.10) 

These rules ensure that signals substantially fill the A/D and D/A converter ranges (thereby 
minimizing quantization noise), while only rarely allowing the signals to exceed those ranges. 

For our system, the converter ranges are Vfs = 1 V. Typical gain values are Kpa = 14.4 V/V, 
and Ka = 0.001 A/V. 



e. Phase-lead Compensation 



The optimal controller described in Eq. (A.l) was derived under the assumption that the re- 



mainder of the control loop in Fig. [s] consists only of the single- mode (2nd order) cantilever model. 
As we have seen, any practical implementation includes additional dynamic components. The pri- 
mary effect of these components is to introduce additional phase lag within the relatively narrow 
bandwidth of the controlled cantilever. Generally, phase lag tends to lower the damping, causing 



the performance to deviate from the predictions of Eqs. (A. 5) and (A. 6) 
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Sources of additional phase lag include: 



• Finite bandwidths of the interferometer, its pre- amplifier, and the actuator coil current 
amplifier. Each device is characterized by low gain and large phase lag at high frequencies. 
Even components with bandwidths exceeding 0.1 MHz, will each contribute a few degrees of 
phase lag at the cantilever frequency. 

• The A/D conversion delay causes an effective phase lag of ioT^i radians. 

• Similarly, the D/A converter zero-order hold, of sample time T, adds ujT/2 radians of phase 
lag. 

While analytical predictions are useful, the most reliable method of determining the total phase 
lag is by careful measurement of the actual system. Minor phase lag will not seriously alter the 
system behavior and, in many cases, can be ignored. Of course, we could have accounted for 
the transfer function of each of these additional components during the design of the controller. 
However, a closed-form solution for the controller would not be possible, and the insight afforded 



by Eq. (A.l) would be lost. Instead, we can easily compensate for the phase lag by augmenting 
the control law with a phase- fead dynamic component. 

Consider the following simple phase-lead transfer function 



Gc{s) 



(1 -I- rjTs) 
+ rs) 



where 



(A.ll) 



V 



1 + sin ( 

i 



and 



sm ( 



As shown in Fig. 10 at the center frequency iVm the magnitude of the transfer function is unity, 
and the phase angle is (pm- By assigning uJm to the controller natural frequency, and concatenating 



this compensator with the controller definition in Eq. (A.l ) phase lags of up to 60 degrees can be 
corrected. 



/. Tustin's Approximation 



Having considered all of the components in the control loop, we derive a practical discrete-time 



equivalent of the control law (A.l). A variety of transformation methods can be used: numerical 
integration, zero-pole matching, hold equivalents, etc. Here we will use a trapezoidal numerical 
integration technique, known in control circles as Tustin's method or the bilinear transformation. 

Given a continuous-time transfer function H{s), the z-transform of a Tustin discrete equivalent 
can be found by substitution 



H(z)=H(s) 



2 z-1 
T z+1 



(A.12) 



where T is the sampling interval. 
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CO 

m 

Frequency (log scale) 
FIG. 10: Phase-lead compensation Bode diagram 

The resulting discrete transfer function H{z) is the ratio of two polynomials in z and is of the 
form 

M 

H{z) = '-f (A.13) 



k=0 

For purposes of controller computation, the constants and bk are arranged into the corre- 



sponding difference equation of the form shown in (A.8). 



g. Prewarping 

The bilinear transformation maps the stable region of the s-plane (the left half plane) exactly 
into the finite region within the unit circle of the z-plane. Clearly, distortion takes place in this 
mapping, such that the frequency responses of H{s) and H(z) differ. We can force the discrete 
equivalent to match H(s) at the controller natural frequency by "prewarping" the continuous-time 
transfer function. Expressing H(s) as H{^), replace lOqc with coa such that 

a;. = |tan'^ (A.14) 

Notice that uja approaches ujoc for high sampling rates. Stated as a frequency substitution the 
prewarped Tustin equivalent is described as 



H{z) = 



Z-1 



(A.15) 



LJoc taXL {id ocT / 2) z+1 
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4. The Complete Controller 



Finally, the complete discrete-time controller can be assembled. To maintain the loop gain 
G{s)H{s) in Fig. |8j we must adjust the original control law (A.l) to account for the additional 
gain elements in Fig. |9j Combining the optimal controller with the phase-lead compensator we 
have the continuous-time control element. 



KaKpaKcKi^ 

(g + Zoc){l + r]Ts) 

(s2 + {uJoc/Qoc)s + a;2J(l + TS) 



(A.16) 



The coefficients of the corresponding difference equation ( |A.8 ) are derived using the prewarped 
Tustin's approximation (A. 15) 



ao 
ai 

as 



i-Lo' - 3r' + 3lo'^ + uj'^t')Qoc + oj 

(w'3 + w'V + a;' + r')goc 
{-oj' + 3r' + 3w'3 - J'^t')Q 



.12 



■ uj'"^ + uj't' 



12 



Jt' 



{J -t' + -U0''^T')Qa 



■ a;'2 + uj't' 



(A.17) 



and 



^0 



bi 



QocJK' 



KaKpaKcKi^ 

uj' + ijt' + uj'rjz'^^T' + J'^z\ 



12 ^1 

oc 



+ (j'V + + t')Qoc + + Jt' 



bo 



UJ 



7]t' + Jr]z'^y + 3J'^z] 



12 J 

oc 



12 r,l 



J -|- r/r' -|- Jtjz'^^t' + uj 
^ J + rjT' + Jrjz'^y - 3J^z '^, 
^ J + i]t' + Jr]z'g^T' + J'^z] 



UJ 



W' + ^'tIz'oc't' 



12 V 

oc 

^ ^oc 



'J + rjT' + Jrjz'g^T' + J'^z' 



where the normalized parameters are: 
K' 



Kqc 



^ z' 



and J 



tan 



(A.18) 



A final note: Occasionally the controller output calculated from the difference equation, 

will exceed the D/A converter range [0, C^a]- Since some converters wrap overranged values to the 
opposite edge of the range, it is essential that the output value be limited: 



y{n)out = max(min(?/(n),Cda),0) 



(A.19) 
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